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Abstract 

We study the elastic deformations that appear due to tidal and cen- 
trifugal forces acting on an elastic sphere in helical motion in a spherically 
symmetric gravitational field, where gravity is considered to be given by 
either a Newtonian or a Schwarzschild background. We review an exis- 
tence/uniqueness theorem based on the implicit function theorem for the 
nonrelativistic case and give explicit solutions to the linearized elastostatic 
equations in both cases. 

1 Introduction 

In this paper, we consider a sphere composed of homogeneous and isotropic 
elastic material in helical motion around a fixed gravitational centre (an elastic 
planet, so to speak, which, like the moon, turns the same side to the gravitational 
center). Gravity is considered as a given background, described either by a 
Newtonian potential or the Schwarzschild solution of general relativity. The 
elastic sphere is therefore exposed to a combination of tidal and centrifugal 
forces, both of which cause elastic deformations. Due to the helical nature of 
the motion it is possible to go to the co-rotating coordinate system, in which 
we face only an elastostatic problem. 

* Supported by Fonds zur Forderung der wissenschaftlichen Forschung in Osterreich, Pro- 
jekt Nr. P20414-N16. 



The organisation of this paper is as follows: in the sections 2-4 we briefly 
review the general theory of relativistic elastodynamics, relativistic elastostat- 
ics and nonrelativistic elastostatics, viewed under one general framework. To 
get reasonable boundary conditions, in elastostatics it turns out it is preferable 
to use the material picture. In section 5, we consider the general principles of 
existence/uniqueness proofs based on the implicit function theorem in elasto- 
statics, as used by Beig and Schmidt in g], 0, [§]. After having considered 
the linearization of the elastostatic equations in section 5, we review, in section 

6, the existence proof given in [5] for the elastic sphere on a circular orbit in a 
Newtonian gravitational field in a version adapted to our purposes. In section 

7, we then give the explicit solution for the linearized elastostatic equations 
in this situation, making use of an ansatz from Papkovich and Neuber ([T3J, 
[12], also see [10] and the references therein). In section 8, we also consider the 
linearized elastostatic equations for the relativistic case and give the solution. 
A existence/uniqueness proof similar to the nonrelativistic one is presumably 
possible, but some details still need to be worked out yet. This will be part of 
a future publication. 



2 Relativistic Elastodynamics 

We start from full relativistic elastodynamics, as laid down in [3J, [TB] or [2J. 
The holds of the theory are mappings 

f-(M, giMV )^(B,V ABC ) (1) 

from the spacetime M, equipped with the spacetime metric (assumed to 
have signature (—,+,+,+)), to the body manifold B. It can be thought of 
as the abstract collection of points ("atoms") the elastic material in question 
is composed of. Generally, B is a open, connected and bounded subset of M? 
with smooth boundary dB, and equipped with a volume form Vabc- Its single 
componend will be denoted V{X) = Vi2s{X) > 0. We will use coordinates X A 
on the body manifold B and x^ on the spacetime M. 

The mapping f A {x^), which is also called configuration, is required to uniquely 
define a future-pointing, timelike, normalized {u^u" g pv — — 1) 4-velocity vector 
it^ defined by 

f A ,^ = (2) 

i.e. the inverse of a point X A g B is a timelike curve in M., the trajectory 
of that point. We further require f A (x^) to be orientation-preserving, which 
allows us to define a positive particle number density n by 

f A A* a )f E 'Ax")f c ', p {x°)V ABC (f D {x°)) = ne^pix^u" (3) 

Using the configuration mappings f A {x^) from M. to B to describe elastic 
materials is called spatial or Eulerian picture. This is the more natural way to 
go in relativistic elastodynamics. It is also possible to use mappings ^(X A ) : 
B — > A4, called deformations, which are in a certain sense inverse to the f A s. 
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This is called material or Lagrangian picture, which is more commonly used 
in non-relativistic elasticity theory, but in full relativistic elasticity theory, the 
spatial picture is much more natural. We will later on introduce and use the 
material picture when dealing with elastostatic situations, though. 

From the configuration gradient f A ^, we can now define the strain tensor 

H AB (x x ) = f A ^(x x )f B ,„(x x )g^(x x ) (4) 

By definition, it is symmetric and positive definite, so it also has an inverse 
Hab- 

An elastic material is specified by giving an stored- energy function 

w = w{f A ,f A ^x») (5) 

which describes the elastic potential energy per particle. We assume w to be 
covariant under spatial diffeomorphisms (the so-called "principle of material 
frame indifference ") , which implies that w depends on the configuration gradient 
f A ltl and x^ just via the strain tensor (see e.g. [3]): 

w = w(H AB J c ) (6) 

From differentiating w, we can define the second Piola-Kirchhoff stress tensor 

duu 

the first Piola-Kirchhoff stress tensor 

v A = f B ^ BC H CA (8) 
and the Cauchy stress tensor 

dw 

The equations of motion for the fields f A are then given as the Euler-Lagrange 
equations of the action 

S= f ne{H AB ,f c )^/- dot g^d 4 x (fO) 
Jm 

where e is the relativistic total energy density given by rest energy plus stored 
energy function: 

e(H AB J c ) = ci+w(H AB J c ) (11) 



<V = n f A ^f B .vOAB = -^n— (9) 



3 Relativistic Elastostatics 
3.1 Material picture 

We now want to consider time-independant situations in a way similar to 0]. 
For this, we suppose that we have an everywhere timelike Killing vector field 
and that the velocity vector field is proportional to it, so that we have 

£(,f A = Cf A ,» = (12) 
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We further suppose that the spacetime M. can be written as product of space 
and time 

M=RxN (13) 

where the space manifold N is the quotient of Ai by the isomctry group gen- 
erated by The metric hij on N is in a suitable coordinate system (t,x l ) on 
IxJV given by 

hij = gij - goigoj/goa (14) 

whereas the inverse metric h lJ is just equal to the spatial part of the inverse 
metric of A4 : 

h ij = g ij (15) 
This is why we immediately get for the strain tensor 

H AB = f A ,J B ^" = f A ,if B ,jhV (16) 

The particle number density n is now determined via the metric epsilon tensor 
tijk on N belonging to hij via 

f A ,i(x l )f B ,j(x l )f c tk (x l )V ABC (f D (x 1 )) = ne m {x l ) (17) 

which is equivalent to 

„" = do, {H DE ) - (18) 
Vabc det hij 



To adapt the action ( 10 1 to the elastostatic setting, we need to split yj— det g^ u . 
To do so, we use the cofactor formula for the inverse matrix, applied to calculate 
goo from g 1 *", which gives 

5oo = ^— det h ij (19) 
det g^ v 



so the action principle ( 10 ) can be rewritten as 



IxN 

we rewrite goo further as 



S= [ n e(H AB J c )y/^m)V det hij <? x dt ( 20 ) 

JRxN 



9ao = c 2 e^ 2 (21) 



we insert this in (20), divide by c and leave the integration over t away to get 
the reduced action 

S= [ ne(H AB J c )e u/c2 ^/feth~ d 3 x (22) 

J N 

By varying this action principle, we get the Euler-Lagrange equations 

e'^Dj (e^W) - n (l + ™) D t U = (23) 
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where Dj is the metric covariant derivative of hij, together with the boundary 
conditions of vanishing surface traction: 

°^ n i|/-i(3B) =° ( 24 ) 

with nij being a co-normal vector to the boundary. Here an awkward feature of 
the spatial picture shows up: the inverse of f A appears in the description of the 
boundary. This difficulty can be avoided by switching to the material picture, 
which we are now going to introduce. 

3.2 Material picture 

In the material or Lagrangian representation of elasticitiy theory, 

: N -> B (25) 
In the elastostatic setting, the <p l are just the inverse of f A : 

^(/V)) = x\ f A (4> l (X B )) = X A (26) 

The map <f> % (X ) is called deformation. From its gradient, we can define another 
form of the strain tensor: 

H AB (X C ) = <f>\ A (X C W, B (X c )%(0(X c )) (27) 

It is the inverse to H AB and as such, contains the same information. 

The first Piola-Kirchhoff stress tensor already introduced in equation {[sj) 
has a particularly simple form in the material picture: 

a i A = f B , i a BC H CA = ^- (28) 

It will also appear in the elastostatic equations and the boundary conditions. 



To get them, we have to transform the action principle ( 22 ) to the form 



S = f e(H AB ,X c )e u / c2 Vd 3 X (29) 
Jb 



From the variation of ( 29 ) , we get the elastostatic equation in the material 
picture: 

e- u ^V A (> C V) - (l + J) diU = (30) 
where is the so-called material derivative given by 



Vao-S = ^d A (Va A ) - a fc V, A (31) 



The r*L denote the Christoffcl symbols of the metric hij, evaluated at cj) m (X B ). 
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The elastostatic equations in the material picture ( 30 ) are accompanied by 
the following boundary conditions: 



where ua is a conormal vector to dB. As already noted, the structure of this 



boundary conditions is much nicer than in the spatial picture (equation (24)) 
because the boundary is fixed and does not depend on the fields. 



4 Nonrelativistic Elastostatics 

The equations governing nonrelativistic elastostatic can be retrieved as limit 



for - — > 0, which can be applied directly to the action principle (29) using the 



decomposition (11) and 

U „ / 1 



e 



W=l+»+Ol±) (33) 



if one neglects a constant "Null Lagrangian" proportional to c 2 , which goes 
to infinity in the nonrelatvistic limit, but doesn't contribute to the resulting 
Euler-Lagrange equations to get 

S= f (w{H AB ,X c ) + U{^(X c ))) Vd 3 X (34) 
Jb 

From varying this action, we get the nonrelativistic elastostatic equations 

V A (Ti A - d t U = (35) 
again together with the boundary conditions 

ai A n A \ gB = (36) 

Alternatively, those could also have been retrieved by applying the limit - — > 



directly to (30), (32) 



5 Existence/Uniqueness Proof 
5.1 Reference State 

For the purpose of finding existence/uniqueness results, we assume that there is 
a stress- free reference state cf> 1 (all quantities referring to this reference configura- 
tion will be denoted with a bar). We will later use the implicit function theorem 
for finding solutions near the reference state. Stressfreeness (i.e. &i A = 0) is 
equivalent to (because of equation 

ID yB^ 

= (37) 

4>=4> 



dw(H CD ,X E ) 
dH AB 
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We can always choose coordinate systems so that the reference deformation 
map is just the identity map on B, i.e. we identify B with the subset of space 
N occupied by the elastic body in the reference state: 



5\X A 



(38) 



In such a coordinate system, H AB = S AB can be considered as metric on B. 
We further assume that the material has a constant density po in the reference 
state, which is equivalent to the volume form on B being the metric one and its 
component V(X) = const. 

5.2 Operator formulation 



We can write equation ( 35 1 in the form 

ei[^'] + W] = o 

where the elasticity operator £i[0 J ] is defined by 

W] = W A a z A 



It is quasi- linear, as can be easily calculated from formula (28): 

d A a A = {AH AE H BF f c 4 f D d L CBDF f c ,iH AB HcD°~j 
-/V/-/ 5 ^) {<P,ab)+1.o. 
where the elasticity tensor Labcd is given by 

d 2 w 

^abcd — qjjabqjjcd 



(39) 
(40) 

(41) 
(42) 



For its value in the reference state, we require the uniform pointwise stability 
condition 

L ABC dM ab M cd > VM AB = M< Afl > ^ (43) 

In case of homogeneous and isotropic materials, the elasticity tensor can be 
expressed using the Lame constants A and \x and the density po of the material 
in the reference state: 



L A 



BCD 



1 

4^ 



(XSabScd + p (Sac^bd + SadSbc)) 



For the Lame constant, the condition (43) is equivalent to 

3A + 2p > 0, fj, > 



(44) 



(45) 



The force term in equation (39) is just given as Nemitskii or composition 
operator: 

T l W]{X A ) = {K l o^)(X A ) (46) 
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where Ki — —diU, as usual. It is also called "load" in the context of elastostat- 
ics. 



It will turn out to be convenient to extend equation ( 39 ) in order to contain 
the boundary conditions as well. We thus introduce operators E,F : 6 — > L, 
which assign tuples (&j,<7j) e £ to deformation maps <p l € 6, where hi is a 
volume force definded on B and <7, is a surface force defined on dB. 

In order for the operators E and i 7, to be well-defined and C 1 , we assume 
the spaces C and L to be the Sobolev spaces standard in elasticity theory (see 
[13] . [Z]) : for C we use a neighbourhood of the reference configuration (j> 1 in 
W 2 ' P (B, R 3 ) with p > 3, which we assume to be small enough so that every 
<fi l € C has an C 1 -inverse. For L, which is also called load space, we use 

£ = W°- p (B, R 3 ) x W^^idB, M 3 ) (47) 

Thus, we now have to solve the equation 

E[<jy>} + tF[(t> j ] = (48) 

where E and F are given by 

EmX B ) = ((V A( x^)(X B ), {o t A n A )(X B )\ dB ) (49) 

and 

FM(l B ) = ((if 1 o^)(l B ),0) (50) 
the second component is just 0, because we always use the boundary condition 



of vanishing surface traction (equation (32)). It is well known that such an 
composition operator is C 1 if the vector field Ki is (see [HJ). The variable t 
will become our linearization parameter. 

An important property of the operator E is that it does not have full range: 
if a load (bi,o~i) G E(G), it has to satisfy the equilibrium conditions for the six 
Euclidean Killing vectors 



(Z°(f>)(X)h{X)dV(X)+ / (Zo4)(X)v i (X)dS(X) = (51) 
B Job 

This is most conveniently seen in the spatial picture by using the symmetry of 
cr^, the Killing equation for £ l and Stokes' theorem. The intuitive meaning of 
this is that the total force and total torque on B have to vanish, otherwise the 
body would be set into motion, and no static solutions would exist. 

Finally let us review some well-known facts about the linearization SE at 
4> l = 4> l (see [H]): it is a Fredholm operator with a kernel given by elements of 
the form 

54>\X A ) = {Co4>){X A ) (52) 

where £ l runs through the six Euclidean Killing vectors. Also the image of 
5E has co-dimension six, and is given by the subspace of L that satisfies the 



equilibrium conditions (51) for 
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We now want to invoke the implicit function theorem to prove existence and 
uniqueness of solutions of ( 48 1 . This is possible, because we have already seen 
that both operators E and F are C 1 in the appropriate sense. We are thus 
considering one-parametric families of solutions (t being the parameter) with 
<f>Q = <j) % . As a first step, we note that we obviously have E[(j) 1 ] — because 
of the stressfreeness of the reference configuration. But unfortunately, since 
the linearization of the operator E at 4> l has a non-trivial kernel and range, 
we cannot apply the implicit function theorem directly. We circumvent this 
difficulty by a two-step procedure: 

The first step is to restrict C and £ to C and L so that the linearization of 
E becomes an isomorphism. To deal with the non-trivial range of E, we define 
£ := E(L) and choose a fixed 6-dimensional complement Sq of it, i.e. 



£ = E(L) © S 6 



(53) 



This defines a unique projection operator P from £ to £ along Sq : 

P : £ © S 6 -> £ (54) 



We now apply P to equation ( 48 1 



'oE\ 



'oF\ 







(55) 



The linearization of (P o E) \<jjj\ is just (P o SE) [Scff], which is onto £ by defini- 
tion. 

To solve the problem of the non-trivial kernel, we consider solutions of the 
form 

4(X A )=Q + 4,i(X A ) (56) 

where £J is an arbitrary element from the kernel of E and c/) 1 is from a comple- 
ment of the kernel of E that contains <jf> 1 , which we will denote C. 

With this restriction, P o SE is now an isomorphism from C to £ for each Q 
hold fixed, and we can apply the implicit function theorem to obtain a unique 
solution for ^(X A ). 

As second step, it remains to show that it is possible to find Q. so that the 
equilibrium conditions (51 1 are satisfied. In other words, one needs to check 



that the full equations (48) rather than just the projected ones (55) are valid 



by a suitable choice of This can then be done using the finite-dimensional 
implicit function theorem. For details, we refer to the concrete example below. 



6 Linearization 

The perturbation of cj) 1 is defined as 



dt 



(57) 



t=o 
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From perturbing ( 16 1 we get the expression 

6H AB = Sf \,f".Ji' ! + f A ,iSf B d h ij = -2Sf8?S(/)^ (58) 
The perturbation for the first Piola-Kirchhoff stress tensor can be calculated 



from (28): 



8ai A = -2H AB f c iL BCDE 5H DE (59) 



Here, it was used that the reference configuration is stress-free by assumption; 
otherwise, additional terms would appear in (59), which would make our treat- 
ment much more complicated. 

Using this, the definition of Labcd (equation (44)) and our assumption on 
the coordinate system (</>*, a — 8 A ), we get 

5<j iA = — (\5 iA dj54P + fx {8<j)i. A + 54>A,ij) (60) 
Po 

which leads to the standard expression of the linearized elasticity operator 

d A 5a A = — ((A + ^didjSclJ + fiA5(j)i) (61) 
Po 

In order to get correct physical dimensions, we define the displacement vector 
as 

u i (x j )=tS4> i (f A (x j )) (62) 
We have thus arrived at the standard equation of linearized elastostatics: 

puij.j + (A + fj)uk,ki + tp K l = (63) 

subject to the boundary conditions 

5o-i A n A \ dB = (64) 



7 Circular orbits in Newtonian Gravity 

We now consider a spherical elastic body with radius a on a circular orbit around 
a fixed center with mass M in the distance L in the case of Newtonian gravity. 
The force density is given by K t = —diU with 

g = _^,-(^H» + ^j (65) 

where m is a unit vector orthogonal to the orbit plane and n is the unit vector 
connecting the body to the force center, and ui is the angular frequency. Ob- 
viously, the force field Ki vanishes on a circle given by r = L if L satisfies the 
relation 

uj L= -jj- (66) 
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On this orbit, we can write the potential as 



y = _^(x-frm)m)' + gj (f , 7) 

We will use w 2 = as our linearization parameter t. 

The whole threatment of this problem can be simplified by utilizing the 
mirror symmetry along the planes orthogonal to n and n x m. In a cartesian 
coordinate system with the axes (— n, n x m, m), this means: 

^(X 1 ,^ 2 ,^ 3 ) = ^> 1 (X 1 ,-X 2 ,X 3 ) = cf> 1 (X 1 ,X 2 ,-X 3 ) 
(f) 2 (X\X 2 ,X 3 ) = -c/) 2 (X\-X 2 ,X 3 ) = c/) 2 (X\X 2 ,-X 3 ) (68) 

3 (a\a 2 ,x 3 ) = 4> 3 {x 1 .-x 2 ,x 3 ) = -(j} 3 (x\x 2 ,-x 3 ) 

We can restrict the spaces C and L of configurations and loads to those also 
satisfying this conditions, and denote them G sym and L S ym)- This is possible, 
because G sym still contains the identity, which we use as </>*. Because of the as- 
sumption of homogeneity and isotropy, the operators E and F can be restricted 

to gO G sym ^ ^sym- 

The only Killing vector that satisfies the conditions ( |68| is the translational 
one n % (pointing from the body to the gravitational centre). This generates 
the Kernel of E : G sym — > £> sym , i-e. it is only one-dimensional. Also, the 



equilibrium conditions (51) are all automatically satisfied except for £ l = n\ 
i.e. also the co-range of E : G sym — > L sym only has dimension 1. We thus have 
to consider just <fi l of the following form 

4(X A ) = C(t)n i + &(X A ) (69) 

However, the first step of the proof from subsection |5.2| still remains valid in 
analogous terms. We thus have existence and uniqueness for (\>\{X ). As second 
step, it remains now to determine C(t) from the equilibrium condition 

N (t, C(t)) = I rj (Ki o ((>) (X)d 3 X (70) 

with Ki being the force field, divided by the linearization constant: 

^=d 4 Q(x 2 -(x,m) 2 )+^) (71) 

We again use the implicit function theorem; but this time, the finite-dimensional 
one is sufficient. First we need to show that iV(0, 0) = 0, i.e. that the force field 
Ki is equilibrated at the reference state <fi l . To do so, we use that the integrand 



in ( 70 1 is a harmonic function 

A (rtKiix)) = n'AK^x) = (72) 
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We further use that B is a ball; we can thus invoke the mean value theorem for 
harmonic functions (see e.g. [5]) to get 



N(0,0)= [ rtKiix^x = \B\ {rtKi) (x ) = 
Jb 

where we have used that the center xq of B, the force field Ki vanishes. 



(73) 



Now it remains to show that f^(0,0) ^ 0. From (70) we calculate (using 
the decomposition (56)) 



dN 
dC 



(0,0) 



. n> ) d 3 X 



(74) 



Again, the integrand is a harmonic function, so the integral can be evaluated 
using the mean value theorem for harmonic functions; an easy calculation shows 

that [{Kin 1 ^ , (xq) — 3, so we get 



dN 
dC 



(0,0) =3\B\ =47ra 3 



(75) 



which is non-vanishing, as required by the implicit function theorem. We have 
thus shown existence and uniquenes of both C(t) and <f>\{X A ), hence <pl(X A ), 
for small values of the parameter t = uj 2 . 



8 Explicit solution 

For convenience, we use the linearized strain tensor, which is given by 



(76) 



We will use spherical coordinates with normalized coordinate basis vectors; in 



such a coordinate system, we get for the components of ( 76 ) 



Hrh, 



du r 
dr 

1 dug 
r~d~6 4 
1 / cos 



r 



H r g = 



h T I 

r 80 dr \ r 



Hg 



sin 

1 (\du. 

2 I 

du<j, 
dr 
du^ 



-Ug + U r + 



1 du^ 
smO d<p 
d_ (Ug 



(77) 



i 

2^ 



!M + _ 

cos 9 



1 duj 



sin i 



] u <t> 



1 dug 

smO dcj) 
1 d 



1 d 

txH = divu = — -— - (r 2 u r ) + . 

r z dr x ' r sm dO 



(ug sin 9) 



1 du^ 
r sin 9 dd> 
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the linearized stress tensor 5<JiA is then connected to the strain tensor via the 
generalized Hooke's law for homogeneous and isotropic materials: 

5a iA = 2fiH tJ S J A + A tr HS iA (78) 

We now use the following strategy to solve the linearized elastostatic equations 



with the force given by the potential (67): we first decompose the potential 



in two parts in subsection 8.1 where both parts are axisymmetric, and one 
being equivalent to the problem of a self-rotating sphere, already solved in [4]. 
Because of linearity, both parts can be threated independantly; we concern 
ourself only with the yet unsolved part of the problem. Using the assumption of 
axisymmetry, it is relatively easy to find a particular solution first (in subsection 
8.2[ ), that is, a solution to the inhomogeneous linearized elastostatic equations 
which doesn't necessarily also satisfy the boundary conditions yet. In a second 



step (subsection 8.3 ), we then find, using again axisymmetry and the ansatz from 
Papkovich and Neuber, the general solution to the homogeneous elastostatic 
equations in terms of a Legendre polynomial series. "General" here means that 
the homogeneous solution contains enough yet unspecified constants so it may 
be matched to any boundary conditions. Then the particular solution and the 
homogeneous solution are added and the constants specified by imposing the 
zero-traction boundary condition. 

8.1 Force decomposition 



The potential ( |67[ ) has the disadvantage of being not rotationally symmetric, 
but it can be decomposed in two parts with that property by a simple coordinate 
transformation: a shift of the coordinate origin to the center of the body at —Ln 
(i.e. by y = x + Ln). The centrifugal term in the potential then becomes 

~\ ( x - ( x > m)m) 2 = -i (x 2 - (x, m) 2 ) = 
2 \ (79) 

= --(y 2 -2L( y ,n) + L 2 -( y ,m) 2 ) 

The term — \ (y 2 — (y, m) 2 ) is just the centrifugal potential caused by a rotation 
of the body along an axis through its center; this problem has already been 
solved (see @], [S])- In the end, this solutions can be just added because of 
the linearity of the linearized elastic equation. The constant term proportional 

j 2 

to =n can be neglected because it doesn't contribute to the force density. The 
remaining term 

L(y,n) = LyP 1 (cos6) (80) 

is the part of the potential that is going to be used further on. It has the 
advantage of being harmonic and symmetric under rotations along the n-axis. 

Since we have y < L, we can expand the Newtonian potential in Legendre 
polynomials: 



1 1 1 ^ 



in| L v /l-2£cos0+£ L n=0 



(!)"p„(cos0) (81) 
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Plugging (81) and (80) into (67), the n — 1 term cancels, and the constant 
n = term can again be neglected. Thus, the effective potential for which the 
static elastic equation is to be solved becomes 



U = -L*l^(£) n p n (cos6)j (82) 

Further on, we are going to rename y to r again for simplicity. 
8.2 Particular solution 

We first search a particular solution (i.e. one that doesn't necessarily fit the 
boundary conditions) for the inhomogeneous elastic equation 

<a + ^r<ki = -rKi (83) 

If we assume that the force is given by a potential (Ki = —dtU), we can find a 
particular solution by also introducing a potential for the displacement vector: 
uf = ^diX- By inserting this ansatz in (83), we get 

^^Ax = d % U (84) 
A* 

which can be integrated once to get 

^=xh u < 85 » 

In the case the potential U is a harmonic function, there can be found a solution 
to this equation easily: we can write it as sum of harmonic polynomials: 

oo 

U(r,e,4>) = Y l U n (86) 

In the case of rotational symmetry along the z-axis, these are proportional to 
the Legendre polynomials: 

U n = E n r n P n (cos 9) (87) 

Since we have 

Ar 2 U n = U n Ar 2 + r 2 AU n + kx&Ui, = 2 (2n + 3) U n (88) 



we can write down a solution to equation ( 85 ) immediately 

oo „ 

X = T ;r77~~ — ~r ^^r n+2 P n (cos 9) (89) 
A ^ 2 (2n + 3) A + 2u v ' y ' 
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We thus get 



^ 2(2n + 3) A + 2// " l 

n— 



E 



i 



2 (2n + 3) A + 2a 

71 — 



r" +1 (-sin6i)P^(cos( 



(90) 



«r = o 



and further (using the abbreviation k n = ^ 2f i 2(2»+3) anc ^ with the Legendre 
polynomials and their derivatives always understood to be evaluated at cos 9) 
for the strain tensor 



Hj. 



Y^k n (n+l)(n + 2)r n P n 

oo 

J2 k n r n (- cos 6P' n + sin 2 6P% + (n + 2)P n ) 

n=0 

oo 

£ fc„r" (cos + (-n 2 + 2)P n ) 

OO 

^ k n r n (- cos ^ + (n + 2)P„) 

oo 

^/c n r"(n + l)(-sin^)P^ 

OO 

trP p = divu p = 2 ( 2 " + 3)k n r n P n 

n=0 

The other components are zero. For the stress tensor we get 



H^ a — 



2(1 ]T k n r n P n ^(n + 1) (n + 2) + - (2n + 3)^ 



n=0 

oo 



n=0 

oo 



A 



2[iJ2 k nr n cos 9P^ + (-n 2 + 2)P n + -{2n + 3)P, 



n=0 

oo 



2/i ^ k n r n - cos + (n + 2)P„ + - (2n + 3)P r , 



2/i53fc„r"(n + l)(-sin0)P; 



For the potential (82 1, we have 



Ei, 



E = E 1 = 
-pLu 2 L 2 ~ n (n > 2) 



(91) 



(92) 



(93) 
(94) 
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8.3 Homogeneous solution 

The homogeneous elastic equation 

/iuf Jj + (A + /i)< w = (95) 
is solved by the ansatz from Papkovich and Neuber 

uf = 2{X + 2fi) B t - 8> (xjBj + B Q ) (96) 
A + fi 

where both Bi and B are assumed to be harmonic, i.e. they satisfy the Laplace 
equation. In spherical coordinates (r, 0, <j>) this becomes 

u? = 2{X + 2 ^ B r -d r (rB r + B ) 

A + fi 

u$ = 2{X + 2 ^ B e -do (rB r + B a ) (97) 
A + fi r 

H _ 2(A + 2^) 
"* A + m * 

Because of the assumption of rotational symmetry along the z-axis, uf and also 
Bi have to be independant of <f>. This is also why vff is just proportional to 
Bq. In cylindrical coordinates (p, 0, z) the components of Bi are more easy to 
handle; the connection to the components in spherical coordinates is given by 

B r = sin 0B p + cos 0B Z (98) 
B B = cos 6B p - sin 6B Z (99) 
5* = B^ (100) 

On the other hand, the connection between the components in cylindrical and 
cartesian coordinates is given by 

B x = cos 4>B p — sin <pB^ (101) 
By = sin <pB p + cos ^B^ (102) 
B z = (103) 



or, in complex notation 



B x +iB y = {B p + iB <t> )e i * (104) 



The components of a harmonic vector are harmonic functions only in cartesian 
coordinates, so B z is a harmonic function, while B p and Bj, are not, but B p e 1 ^ 
and B^e 1 ^ are, so they have to be proportional to the spherical harmonics 
with m = 1, and so B p and B$ arc proportional to the associated Legendre 
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polynomials with m = 1. Thus, they can be expressed by 

oo 

B p = A n r n (- sin 6)P' n (cos 9) 

n=0 
oo 

fl* = £C£r n (- S in0)i*( C OS0) 

OO 

B z = ^A*r>\P n (cos0) 



(105) 
(106) 
(107) 



inserting (1051 into (98) yields 



B r = J2 -A n r n P^(cos 9) sin 2 9 + A* n r n P n (cos 9) cos 9 (108) 



n=0 



P e = ^ -r'M n P;(cos 0) sin 9 cos 9- A* n r n P n (cos 9) sin (9 (109) 



n=0 



using the relations (see e.g. [T7] ) 

(x 2 - 1) P' n {x) = nxP n (x) - nP„_i(x) 

this becomes 

oo 

B,. = A* cos 6> + Y r " (-A„nP„_i + cos 9P n (A n n + A* n )) 

n=l 
oo 

B e = -A* sm9 + ~r n sin 6 (AnP^ + P n (A n n + A*)) 

n=l 



(110) 

(111) 



(112) 
(113) 



One of the four harmonic functions in this ansatz can be chosen arbitrarily, so 
we can significantly simplify everything by setting 



A n n + A* = 



(114) 



using this, performing an index-shift and renaming the A n+ i back to A ni we 
get 



B r = ^ -A n (n + 1 )r n+1 P n (cos ( 

oo 

Bg = J2 ~A n r n+1 sin9P^{cos9] 

n=0 

for the fourth harmonic function Bq we use 

oo 

B O = J2- B nr n Pn(cOS0) 



(115) 
(116) 

(117) 



n=0 
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plugging all this into the ansatz (97), we get 

OO / 

E [Mn + l) 



,H 



,H 



2/i 
A + /i 

3A + 5^\ +1 



n+l 



Bnnr 11 - 1 P„(cos 



53 (A n [n+^^^y^ 1 +B n r n - 1 ) (- sin )./>',( cos 0) (US) 

OO 

53 C„r"(- sin (cos 0) 



Tt = 

this gives the strain tensor 



n=0 

OO 



• / 2^ 



^ = E ( - ( 4. K + + xtm) (» + !) ^ + B ^ 

3X + b ^r n + B„r n -A cos 9PL 



n=0 



At. n 



P, 



n=0 



+ ( A n (n + 



A 


+ i" / 




2/i N 


, — 


A + ^/ 


3A 


+ 5^\ 


A 


+ /•* / 







r n + B n nr n - 2 j P n + 
r n + B n r n - 2 ) (-cos0)P; 



n=0 



oo 

n=0 
oo , 

H H = E C « rn_1 COS0P 'n - + 



trP ff = divu H = E 



n=0 



2/z 

A + fJL 



A n (2n + 3){n + l)r n P n 



(119) 
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which in turn gives the stress tensor 
Sag. = 2/, J2 U n (n + l)(n 2 -n- ) r" + B n n(n - l)r n ~ 1 

+ (A n (n+ 3X + 5 ^~« 



^r n + B n r n - 2 J cos 9P' n 



&4 = 2 M £( ( ^ 2 + 2n - ) r w + B„(n - l)r"- 2 ) (- sin 6)P' n 



n=0 

OO 



oo , 



n=0 



(120) 



To get the overall solution, one has to add the particular solution (90) and 
the homogeneous solution (1181 and submit it to the boundary condition to 
determine the constants A n , B n and C n : 

= 5a rr \ r = a + Sa^ r \ r=a 

{) = 8a%, =a +8a^\ r=a (121) 

= Safyr^a 

The equation for Sa r( f, becomes 

(n - 1)C„ = (122) 

so the C n for n ^ vanish, while C\ is arbitrary. The corresponding displace- 
ment 

W = -<7irsin0 (123) 

corresponds to rigid rotations along the z-axis and can therefore safely be ne- 
glected. 

The boundary conditions for Sa rr and 8a r e become a linear 2x2 system for 
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the variables A n and B„ 



n z + 2n — I a n A n + (n- l)a n - 2 B n = -k n a n (n + 1) (124) 

A + /V 



(n + 1) [ n 2 - n - 3A + 2 ^ ) a ™ A n + n(n - l)a n - 2 B n = (125) 

= -fc n a" f(n + l)(n + 2) + ^(2n + 3)^j (126) 
with the determinant 

D n = (n 2n(n - 1) + 3A + 2 ^ (2n + 1)^1 a 2 ' 1 " 2 (127) 

V a + m / 

Because of > 0, we have D Q < 0, D 1 = and D„ > for n > 2. 

Because fco = k\ = 0, this gives Ao = Bq = A\ = 0, while B\ is arbitrary. The 
corresponding displacements 

u r = Bi cos 9 

R • fl (128) 

ug = —B\ sma 

are rigid translations along the z-axis, lying in the kernel of the linearized elas- 
ticity operator. Thus, the value of B\ has to be calculated from the equilibrium 
conditions after the rest of the solution has been determined. 

For n > 2, we get the following unique solution for the system (124 1: 



A n = -^-k n a 2n - 2 {n - 1) (2(71 + 1) + -(2n + 3) 

1 «. 2 n A + 2/i , N (,_ , ^ , A, 



B„ = -— fc„<T% -n(2n + 3) ( n + 1) + _( n + 2) 

combining everything gives the solution of the elastic equations 
/ / r \ 11+1 a 2 /r\ n ~ 1 \ 

U r = sJ2(Fn(j) +Gn T ^( T ) )P n (cO S 0) 



(129) 



n=2 
00 



n=2 

= 



"+ 1 r a 2 /r\«-!\ dP n , „. (130) 
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with 



pGM 



F n = - 



2 (A + n) 2{X + p) 
n (n + £(n + l)) 



2 „(„_i) + 3^(2n+l) 



G n — nl n 
H n - 



n + 2 + ±{n + 3) 



(131) 



2n(n-l) + ^(2n+l) 

7l(Wl+ |(7l + 2)) 

(n-1) (2n(n-l) + ^(2n+l)) 



It is noteworthy that the Lame constants only occur in ratios except in the scale 
factor s. 



The overall solution ( 130 ) has the strain tensor 



+ (n - 1)G„ 



n=2 



a 

J? \L 

2 



Pn (COS ( 



Hgg = -J2 [(Fn - n(n + l)H n ) i^-j -„ 2 / n -^(-J JP„(cos 



Ha 



n=2 



H„ 



L L 2 \L 
G n 

L 2 



n-2 



cos 9P' n (cos ( 



L 2 VL 



n-2 



P„(cos#) + 
(£f 2 ) (-cos0)P;(cos0) 



rO 



trH 



^(i(F„+n J ff„)(^)' i + (n-l)/ n | 

oo 



f \ n—2 

L 2 \L 



sm9)P' n 



n—2 



(132) 



From that it is easy to check that the boundary conditions (121 1 are indeed 
satisfied using the algebraic identities 

F n +nH n + 2(n-l)I n = Q (133) 
2/i ((« + l)F n + (n- l)G n ) + A ((n + 3)F„ - 3n(n + l)H n ) = (134) 
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Also the elastic equations can be checked: 



(AS) r = 



r 2 dr 



, du r 
dr 



1 8 2 



r 2 &e 2 



(Au) 6 



s 

U 2 

n=2 

r 2 dr 



1 cos 9 du r 
r 2 smO d9 
n-i 



2u r 



d(ug sin 0) 



^ (2nF n + 2n(n + !)#„) P n ( 



, ^ 1 8 2 U0 ^ 1 cos9du 



dr J 



u 



r 2 sinO d6 



r 2 sin 2 i 



2 chx r 

7 2 ~m 



s 

L 2 



J2(2F n + 2(n+l)H n ) 



dd 



(cos 9) 



(AS), = 



(grad div u) 



s 

L 2 



(n(n + 3)P„ ~n 2 (n + l)H n ) (-) P„(cos 



n=2 

oo 



(grad div u) B = ^ ((n + 3)F„ - n(n + !)#„) 



n=2 



dP„ 
~d9 



(cos 9) 



(135) 



(136) 



(grad div u), — 
From the algebraic identity 

fi (2F n + 2(n + l)H n ) + (A + t i) ((n + 3)F n - n(n + l)H n ) 



-2(A + /i) (137) 



it then follows that the solution ( 130[ ) indeed satisfies the elastic equations with 
the potential (82). 

8.4 Displacement 

To complete our solution, we still have to calculate the displacement factor C(t) 
of the Killing part of the solution (compare equation (69)). We do this in first 
order in the linearization parameter t = uj 2 as well, using the explicit formula 
the implicit function theorem provides us with on the equilibrium condition 

dC 



(70): 



dt 



t = - 



1 (Q,oy 2 
f§(o,o) 



(138) 



The denominator has already been calculated in equation (75). It remains to 
compute 



dN 
~dt 



(0,0)* 



A',;/') . u 3 ) d 3 X 



with Ki being the force field, divided by the linearization constant: 

L 3 ' 



Ki = & 



^(x 2 -(x, 



(139) 



(140) 
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The centrifugal term in (i^n 1 ) is just proportional to rij = Pi(cos^)eJ + 

P^ (cos 0)e 8 j, so it doesn't contribute to the integral because of the orthogonality 
conditions for the (associated) Legendre polynomials (see equation (148) below), 
since the sums for n? only start with 1 — 2. So it remains to calculate the 
contribution from 

£3 =L 2 Y (^) Ptjcose) (141) 

From differentiating with resprect to y l , we get 

7-3 QQ / i I \ (-1 

di j, , 2 =. , = T2 =lJ2(-t) (^(«Mff)ef + (-sin0)P/(cos0)e? 
^|y| 2 - 2£.|y|cos0 + £ 2 T^\\ L ) 

. (142) 

The translational Killing vector in direction of the gravitational center is 

C = e\ = cos 0e* - sin #4 (143) 



Inserting this in ( 142 1 yields 



g« 7 , 2 „ ', „ r2 r = (f) (icos0P,(cos0) + sin 2 0P/(cos0)) 

V lyl - 2L|y|cos6> + L 2 J^l 

(144) 

Now using the Legendre polynomial relation 

{x 2 -l)Pl(x)=lxPi(x)-lP l _ 1 (145) 
this becomes (after an index-shift) 

j 3 00 / 

di /i 12 on i irm e = £ E( / + 1 )(f) p '( cos0 ) ( 146 ) 

Vlyl - 2L|y|cos(9 + J L 2 ^ W 

which becomes, after a further derivative 

/ r3 \ 00 / Z— i 



/ r \ dP \ 



(147) 

This can be easily integrated over B, using the following orthogonality conditions 
for the (associated) Legendre polynomials (see [T7]) 

2 

Pi(x)P n (x)dx = - — —Sin 

271 + 1 (148) 
P?\x)p£\x)dx= 2 ^±^ 5ln 
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the integral (1391 then becomes 



dN 
~dt 



(0,0)f = 4ira 3 sY^ 





i 2 ™ n(n - 




(I 


' 2n-\ 





2n 



(149) 



Thus, for small values of the linearization parameter t — uj 2 , we get the 
linear approximation 



C 



dC 



di 



(0,0) W 2 



t=0 

oc 



9JV r(0,0) 



dC 1 



n(n + 1) / F„ + (n + l)i?„ 



2n + l 



2n + 3 



(150) 



Using the formulas for .F n , i?„ and I n , it is easy to see that each summand in 
( 150 ) is positive, i.e. the overall minus sign makes C negative. 
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9 Relativist ic case 

We now want to consider an elastic sphere on a circular orbit in the Schwarz- 
schild metric 

9fl „dx»dx v = - (l - ^) c 2 dt 2 + (l - ' dr 2 + r 2 (d9 2 + sin 2 9d V 2 ) 

(151) 

We first recapitulate that these orbits (for point particles) can be found using 
the following neat property: we assume £ M to be a timelike Killing vector field, 
i.e. it satisfies the Killing equation 

V^ + V^ = (152) 

Then a orbit 7 of the flow of £ M with the normalized tangent vector (4- velocity) 

un = , h (153) 
" V-^p9 Xp 

is geodesic, i.e. satisfies 

u^V^u v = (154) 
if and only if the gradient of £ 2 vanishes everywhere on 7: 

V»(U P 9 Xp ) = 2£ A V„6 = (155) 



This can be easily seen by inserting |153|) into ( 154 1 to get 

L ' 1 r v M ^= = -lr^y,e + r^eev^x (ise) 



and noting that the second term vanishes identically because of the Killing 
equation (152 1. 

We now apply this to the helical Killing vector (remember that constant 
linear combinations of Killing vectors are Killing vectors themselves) 

ed» = d t +ud (p (157) 



of the Schwarzschild metric (|15l|). For its normal value, we get 



< \P 9Xp = I 1 - ^) + u 2 r 2 sin 2 9 (158) 



From forming the gradient of (158), we get the equations 



2(UI ±j 2 r>h> 2 (159) 



2w 2 r 2 sin 9 cos 9 = u 2 r 2 sin 29 = (160) 
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Equation (1601 has the solutions 9 = 0, 9 = ir/2 and 9 — n, but equation 
(159) can only be solved for 9 = ir/2. We thus got the circular geodesies in the 
equatorial plane with r = L = const, where tu, M and L have to satisfy the 
same relation as in the Newtonian case: 

u 2 L = -p- (161) 



In order for £ M to be timelike, i.e. equation ( 158 ) to be negative, it is necessary 
that L > 3GM/c 2 . 

We now consider an elastic sphere moving on such an orbit. We do so by 
changin to the co-rotating coordinate system. We do this by introducing co- 
rotating coordinates: 

ip' = ip + Lot (162) 

Using this coordinate transformation (and renaming ip' back to ip afterwards), 
( 151 ) becomes 

, „ 7 „ ( 2GM oj 2 r 2 sm 2 9\ 2l2 / 2GM\ _1 l2 
g^dx^dx" = - ( 1 - - 2 J c 2 dt 2 + ( 1 - j dr 2 + 

+ 2wr 2 sin 2 9d(pdt + r 2 (d9 2 + sin 2 9 dip 2 ) 
We introduce the potential term 



(163) 



2u 2GM w 2 r 2 sin 2 < 



^=!-^r — ( 164 ) 

The metric hij on N, the quotient of the Schwarzschild spacetime M along the 
helical Killing vector £ M given by (1571, is then given by (see [4]) 

hij = gij - goi9oj/goo (165) 



With (163), we get 



2GM\ 



-l 



//,: ,-,/.,•'./.,-' ( I ) flr'~' f r 2 d9 2 + ( r 2 sin 2 6 + ^r 4 sin 4 0e~^ ) dip 2 

(166) 



\ c z r J \ c z 

The elastostatic equations are 

e^Vife?^) - (l + ^)diU = (167) 

with 

V A a A = V~ 1 8a (Vai A ) - T k ij( j k A & A (168) 
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The Christoffcl symbols of the metric ( 166 ) become 



2 1 c 2 r 2 \ c 2 r 

1, 2uj 2 L 3 

t 86 — —~nge r /h, rr — — r H ^ — 

2 c z 

r (^j^j — ~ hip(p,r I hyp t sin ~\~ CD 

T e re = \hee,r/h B e = - (169) 
2 r 



r w = —~h Wi o/hee — ~ sin#cos# + 



,2 



r v e v = ^h^g/htptp — cot9 + o 



The not mentioned ones are zero. For the last two Christoffel symbols we have 
used the Taylor expansion 

^ = ^ = ^( 1 -^ si " 2 ^ + °(^)) (170) 

We notice that the Christoffel symbols consist of two parts: the regular Christof- 
fel symbols of flat space in spherical coordinates (which we will denote T l jk from 
now), plus some correction terms proportional to u; 2 /c 2 , which we will denote 

Like in 0], we treat the relativistic equation by splitting it in the nonrela- 
tivistic problem plus some correction terms. Thus we write 

e ' £Tj = £Tj + U) (Ji [J. I 1 J 

Inserting this decomposition and (1681 in (167) and multiplying with e u ^ c 



yields 



d t U = 
(172) 

The linearized equation can then be obtained by differentiating with respect to 
w 2 and setting it zero afterwards. For simplicity, we perturb around an relaxed 
configuration, i.e. we assume 

<vV=m=*)=° (173) 

As usual, we will identify B with the part of physical space N occupied by the 
body in the reference configuration, i.e. 

/V) = Sfx* (174) 
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The superscript ?? will be used from now on to denote quantities referring to 
the reference configuration only. We also assume the stored energy function w 
to be zero in the reference configuration 



w(H AB ,X)\ {u2 ^ $) =0 (175) 

this is necessary, because w explicitely occurs in the relativistic elastostatic equa- 
tions (as opposed to the nonrelativistic ones), and thus they are not invariant 
under adding a constant to w. 

The perturbation of &i A is then just equal to the standard expression well- 
known from nonrelativistic elasticity theory: 

Sbi A = -2H AB f c i L BCDE SH DE (176) 

where the Lbcde ar e given by the Lame constants: 

Lbcde = j (XSbc?>de + V {Seduce + SbeScd)) (177) 

The force term just becomes the nonrelativistic one in the perturbation process 
(remember that w — and e u ^ c — 1 in the reference configuration): 



d 



d{u 2 ) 



K"(i + 5)«0 



L '"' ' V i -:/. <1T« 
" 2 =o ' ' V o 



For the Christoffel symbol term in (172 1, we get 
d 



d{uj 2 ) 



(~e u l c2 T\^ A a k A ) 



T k l3 5\{5a A +a A l 2 ^) (179) 



The Christoffel symbols f k ij together with the partial derivatives form the reg- 
ular affine connection on flat space. When we denote this Va, we get for the 
perturbed elastostatic equation 

V^i Jt + V A of i y *| ua=0 + / i = (180) 

The first term is the nonrelativistic elasticity operator, which becomes with 
( 176] ) and ( fl77j ) 

V A Sa A = fiASft + (A + //)graddiv<^ (181) 
The perturbations 8(f) 1 are proportional to the displacement vector u l : 

u\x j )=u 2 8(j> i {f A {x i )) (182) 
The perturbations of the correction term of the stress tensor is given by 

a A \ u2 _ = lim \ (a k A - a k A ) = -2H AB f c \L BCDE k DE (183) 
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all other terms vanish because of the assumption of stressfreeness of the reference 

(184) 



configuration (173 1, and K DE is defined by 

d 



K 



AB 



-H 



AB 



the inverse of the three-metric (|166|) inserted into (|184|) gives 





2L J 

c 2 r 



o o-i 



(185) 



and the same components for K AB because of f A i = Sf. With this and (177), 
equation (183 1 becomes 



1 



u 2 =0 



Atr K5f - iiK AB 5 m 



with the divergence 



(186) 



(187) 



For further use, it is convenient to write JsT AB in cartesian coordinates, which 
we are free to do, since the linearized problem is defined on Euclidean space. 
Then (185) becomes 



k AB = - 



2L 3 x A x B 1 



d 



d<f) J \dcj) 



d 



with the rotational Killing vector 

d_ 

dip 



r ABC 



m B X c 



(188) 



(189) 



where m is the unit vector orthogonal to the rotation plane. 
9.1 Solution 

We now proceed as in the nonrelativistic case by shifting the coordinate origin 
to the center of the body at —in (n is the unit vector pointing from the body to 
the gravitational centre, i.e. (—1,0,0) in the cartesian coordinate system used 
above) by 

y = x + £n (190) 

Thus, y is the normal vector to dB = {|y| = a}, so the boundary conditions to 
the linearized elastostatic equation ( |180 ) to be solved become 



(<5<7i 



A , ~ A\ 



ISA 



|y|=« 



(191) 
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This system is best to handle if one treats the parts of K independently by 



gj^AB + rfi-AB 
X A X l 



AB 



1 







(192) 
(193) 

(194) 



Also, tTj" 4 1 ,,,2 _ n is decomposed in a similar way by inserting 9K AB and r K AB 
in equation (1861. Thus, the elastostatic equation (180) with the boundary 
condition (191) decomposes into: 

1. The non-relativistic equation with the homogeneous boundary conditions: 

V A 5a A + fi = (195) 



S& 



aVa 



|y|=o 



= (196) 



This has already been solved 

2. a relativistic correction term for 9 Oi A \ 
part of the curved metric (166): 



, coming from the gravitational 



V A 5a A +W A g a A \ 



(Sa. 



A _j_ g~A\ 

u 7. 



; 2 =0 



VA 



= 



|y|=» 



(197) 
(198) 



and 



3. a relativistic correction term for r ai A \ 



of the same form as (197), 



(198), coming from the rotational part of the metric (166). 
9.2 Gravitational part 

We consider the relativistic corrections due to the term 

2L 3 X A X B 



Its trace is given by 



so we get 



9R AB = -'- 



tr 9 K 



ga i A = ^ 1 S A 
cr r 



2L 3 



2/iL 3 XiX 1 



Cr r J 

The divergence of 9 Ui A leads to the correction term of the force 

L 3 (2fi - A) x 

t : : - \ - I ; 

r 



g-xA\ 



(199) 

(200) 
(201) 

(202) 
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which can be expressed as the negative gradient of the potential 

U = L ' 3(2M 2 ' A) - = f>|yr^costf (203) 
c r ' — ' 

n=0 

with the coefficients 

E n = \(2^-\)L 2 - n (204) 

With these, one gets a particular solution that is of the same form as in the 
nonrelativistic case: 

P £ 3 2 M -A^ n + 2 f\y\\ n+1 D , m 

M ^^27rTA^2(^T3)Uj Pn(cos0) 

71=0 V / \ / 

P l 3 — a ^ i nyiy +1 , . m ( 205 ) 



c 2 2^ + A ^ 2(2n + 3) V X 

< = o 

Let us us now consider the boundary conditions. The unit normal vector to 
B is yA/a. Inserting this into 9 Oi A yields: 

:= y A a = — ! 5 (206) 

c z r a c z r 6 a 

In order to plug this in our ansatz, we need to develop this in suitable (associ- 
ated) Legendre functions of the form 

oo 

9 °\\=a = £M rl (|) n p„(cos0) (207) 

OO 

9 ^l|y|=a = J2 N n(lY P n(oo S 0)(-sme) (208) 

n=0 

The most convenient way to accomplish this is to start by the generating func- 
tion 

i 1 1 00 / 1 I \ ™ 

- = ; = T y(!v P n (cos9) (209) 

r v /|y| 2 -2i|y|cos^ + L2 L^\L) 

and differentiate it with respect to |y|, which yields 

Lcostf-|y| 1 ^ (W^' ' 



J2 n (-Ti p n(cos9) (210) 



v /|y| 2 -2L|y|cos0 + L2 3 ^ 
The components of 9 Ui in spherical coordinates then become (again, we rename 
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|y| to r from now on to save some writing effort) 



XL 2 c " 



•a 



r \\y\=a 



^(-i P,lr<>„/ 



2^L 2 fa 



71=1 

-2 00 



(211) 



|y|=o 



a * e «l|y|=a 



= 



sinfl^n^ — J P n (cos( 



To get 9 a r in the desired form (2071, we use the identity (see [17] 

1 



xP n (x) 



2n+ 1 

which holds for all n > 1, to get 



((n + l)P n+ i+nP„_i) 



(212) 



2^L 2 



l|y|=« 



E 



V.n=0 



A 

27* 



(cos 6») 



; P„ +1 (cos0) + 



n=l 



2n + l 



2n+ 1 



P„_i(cos( 



(213) 



After two index-shifts, it becomes into the form (|207| with the coefficients 
M, 



2[iL 2 ( A /£\ n(n-l) (n + 1) 2 



2n- 1 



2n + 3 



(214) 



In an analogue way, for 9 ag, we use the identity (also for all n > 1 and taken 
from [T7] as well) 

Pn{x) = {P' n+1 P;_0 (215) 

to get (after two index-shifts as well) it into the form ( |208[ ) with 



2/i£ 2 



L\ n-1 n+1 
a / 2n - 1 2n + 3 



(216) 



We now have everything ready to plug into the ansatz from Papkowitsch 
and Neuber, similar to the nonrelativistic case. The particular solution is given 
by the same formula as in the nonrelativistic case, but with the coefficients 



(204). Also the homogeneous solution is the same as in the nonrelativistic case. 



The free constants A n , B n and C n have to be determined from the boundary 
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conditions 



= { a r8 + °re 

0= K)U 



u / I r— a 



(217) 
(218) 
(219) 



From the boundary condition for and the discrete symmetries imposed, we 
conclude that = 0, as in the nonrelativistic case. The required stress tensor 
components for the particular solution are 



2/i 



L 2 2/i - A 
^2fi + A 



^ f(n + l)(n + 2) A \ /r\« _ . 
g 2(2n + 3) + 2» ^ 



L 2 2/i-A^ 



c 2 2/i + A <^ 2(2n + 3) 

n— 



2(2n + 3) ' 2/i y 
1 / r \ ™ 

--( z ) (-sin0)i*(cosf 



(220) 
(221) 



With this and (|20JJ), (|208j, (|214j) , ([216]), the other two boundary conditions 
become 



(n + 1) n - n - 



3A + 2/i 



A + /t 

L 2 /a\« /3A + 2/i n 2 -n-4 
Z/ I A + 2/i 2(2n + 3) 





2 n(n — 






' 2n- 





(222) 



n 2 + 2n - — ^— ) a"A„ + (n - l)a n ~ 2 5„ = 

L 2 / a\ n f 3\ + 2/i n + 1 f L\ 2 (n 



LJ \ A + 2/i 2(2n + 3) V a / 2rl 



(223) 



For n = 0, we actually have just one equation, because the are component of 
both the particular and the homogeneous solution vanish identically, and Bq 
does not occur in the overall solution, we thus get 



(224) 

3 <r A + 2/i 



which leads to 



< = ~\Kt (225) 
3 cr 

= (226) 

This means that we get a contraction of the elastic sphere due to relativistic 
effects. 
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To get a solution of the system (222), (223) for n > 0, we view the determi- 

(227) 



nant of the system. Like in the non-relativistic case, it is 

D n = -(n -1)( 2n(n - 1) + ^±^(2n + 1) ) a 2 "" 2 



X + fi 



D\ = 0, so there is no unique solution for n = 1, but we see that equation (222 ) 
becomes —2 times equation (223). Thus there are solutions, which are given by 

1 L A + n 



5 c 2 A + 2[i 



(228) 



while i?i is arbitrary. After adding the n — 1 term of the particular solution 

r r 2 cos# + f?i cos( 



( 205 ) to the homogeneous solution, we thus get 

1 L 



10 c 2 
7 L 

10 c 2 ' 



2 (-sin 6>) + Bi(-sin6») 



(229) 
(230) 



The term proportional to B\ is the translational Killing vector e*, which has 
to be determined from the equilibrium conditions after the rest of the solution, 
just as in the nonrelativistic case. 

For n > 1, we have D n < 0, so we can solve the system (222), (223) directly 
to get 



A n 



1 



\+2fj, 



B n = - 



2 L ' n+1 (2n + 3) (2n(n - 1) + ^±^(2n + 1)) 
L 3 1 



c 2 L"- 1 



1 

2n- 1 



2 (2n(n-l) + ^(2n + l) 



(231) 
(232) 



After adding the particular solution (90) and the homogeneous solution (118) 
using the constants (231) and (204), we thus get as overall solution 

, n+l n 1 



n=2 
oo 



END + G 4 



U0 



-!\ dPn 
dd 



(COS 6>) 



(233) 



U0 = 



35 



with 



ur 



'L :i 



F, 



(n + 2) (2n 2 + 



A _ _ 3A+2p 
A+jU 2(A+ M ) 



(2n + 3) ( 2 n(n -l) + ^(2n + l)) 



H n = 



2n 2 + m±8M n+ 1K3A+2M) 
A+^i 2(A+/j) 

(2n + 3) ( 2 n(n-l) + ^±M(2n + l)) 



L 



A+M 



a 7 2n — 1 



(2n(n-l) + ^gi(2n + l) 



(234) 



9.3 Rotational part 

We consider the relativistic corrections due to the term 



r K AB = -^(m x x) A (m x x) 



1 



(235) 



Its trace is given by 



tr r if = -^(x 2 -(x, m) 2 



(236) 



For the sake of completeness, we collect the terms ( |186[ ), (235) and ( 236[ ) to get 



r~ A I 

°i ,„2-_ 



A ±(x 2 -(x,m) 2 ) 



(237) 

To calculate the force correction term in the elastostatic equation (180), we 
calculate 



1 



1 



X[tv r K\ = A-j- (xi — (x, m)mi) 



and the divergence of r K is given by 



1 



V B r i^ = - 2^-(m,x)m y 



Thus the correction term in ( 180 ) becomes 

V7 r~A\ X — n 

Va L»=o ~ 



- (x, m)m 4 ) 



After the coordinate origin shift (190), we get 

1 



r K AB = - 4r (m X (y - Ln)) A (m x (y - Ln)) J 



(238) 



(239) 



(240) 



(241) 



36 



This can be even further decomposed into 

rl K AB = -\(mxy) A (mxy) B (242) 
c 2 

r2 K AB = -4 (-(m x n) A (m x y) B - (m x y) A (m x n) B (243) 

c 

+L(m x n) A (m x n) s ) (244) 

The rl K AB part leads to ( in a cartesian coordinate system with m as the z-axis) 
the system 



//A<i<>' + (A + // igvaddiv<v/ + ( |0+ ) | ,r j = 



W/im + ^ (y 2 - (y, m) 2 ) y 




(245) 



|y|=a 



This describes spheres in rigid rotations along an axis through their center; a 
solution is given in [3]. 

It remains to find a solution to r2 K . By looking at it the right way, 
it becomes immediately clear that it can be considered as the strain tensor 
belonging to the followin perturbation: 

5f A = — ^ ( — (m x n) A (m x n) 3 yj — (m x n.yyj(m x y) A ) (246) 



i.e. 

r2 K AB = <\f A , !' : y)" + f A l 5f B 3 8 13 (247) 

as can be easily seen by elementary vector calculus (Remember that the flat 
metric <5 y is used instead of the curved one , because after linearization, the 
whole problem is considered to be defined on Euclidean space). Thus, the dis- 
placement belonging to minus the one given in equation ( 246 1 has a strain tensor 
that cancels out r2 K AB 1 and thus automatically satisfies both the elastostatic 
equation and the boundary conditions belonging to r2 ai A . So the displacement 
vector is basically identical to (246) (going from Sf A to 8(f) 1 causes another 
minus sign) apart from the factor uj 2 : 

^ ^ (m x n) 1 (m x n)^ yj — (m x nYyj(m x y) 4 ^ (248) 

In a cartesian coordinate system with the axes (m, n x m, n), this means: 

W) = ^y{ ) (249) 

Note that this term breaks the rotational symmetry. 
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Figure 2: A graphic display of the solution including all relativistic correction 
terms, viewed in direction m. The material and orbit parameters are the same 
as for the nonrelativistic picture [l] on page 25 The material used is steel (A = 
108 GPa, (j, = 78 GPa, p = 7860 kg/m 3 ). For the orbit velocity luL = 9.3 km/s 
was used, and L — 3a. The gravitational centre is in the direction of the top 
side of the picture. In order to see the relativistic effects well, the speed of light 
was set to c=20.8 km/s. 
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